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ABSTRACT 

We consider the combined role of the thermal ionization, magneto-rotational and gravitational insta- 
bilities in thin accretion disks around supermassive black holes. We find that in the portions of the disk 
unstable to the ionization instability, the gas remains well coupled to the magnetic field even on the cold, 
neutral branch of the thermal limit cycle. This suggests that the ionization instability is not a significant 
source of large amplitude time-dependent accretion in AGN. We also argue that, for M 10~ 2 M yr -1 , 
the gravitationally unstable and magneto-rotationally unstable regions of the accretion disk overlap; for 
lower accretion rates they may not. Some low-luminosity AGN, e.g. NGC 4258, may thus be in a 
transient phase in which mass is building up in a non-accreting gravitationally and magneto-rotationally 
stable "dead zone." We comment on possible implications of these findings. 

subject headings: accretion, accretion disks - black hole physics - MHD - turbulence - galaxies: nuclei - quasars: 
general 



1. INTRODUCTION 

Thin accretion disks around compact objects are locally 

' unstable at central temperatures ~ 10 4 K, due to large 
opacity changes when hydrogen is partially ionized (Meyer 
& Meyer-Hofmeister 1981). It is generally believed that 
this instability drives the outbursts of dwarf novae (DN; 
Smak 1984; Cannizzo 1993 and references therein) and soft 
X-ray transients (SXTs; van Paradijs & Verbunt 1984; Mi- 
neshige & Wheeler 1989; Dubus et al. 2000). Physically, 

• the disk cannot accrete at a steady rate because a range 
of radii are thermally and viscously unstable. As a result 

, the disk alternates between two stable branches (see §2): 
(1) a hot, ionized branch corresponding to large accretion 

\ rates (outburst phase) and (2) a cold, neutral branch cor- 
responding to low accretion rates (quiescent phase). 

Global time-dependent calculations of the thermal limit 
cycle show that, in order to reproduce observed DN out- 
bursts, viscosity on the cold, neutral, branch must be 
much smaller than on the hot, ionized, branch (a co id ~ 
0.01 << a ho t ~ 0.1; e.g., Smak 1984). For a co i d w a hot , 
the disk only undergoes small amplitude luminosity fluc- 
tuations, rather than the large amplitude outbursts ob- 
served. Gammie & Menou (1998; hereafter GM) showed 
that the magnetic Reynolds number on the cold branch 
in DN disks is Rem w 3 x 10 3 . Numerical simulations of 
the magneto-rotational instability (Balbus & Hawley 1991; 
1998) show that, at such low magnetic Reynolds numbers, 
MHD turbulence and its associated angular momentum 
transport are significantly reduced (Hawley, Gammie, & 
Balbus 1996; Fleming, Stone, & Hawley 2000). Menou 
(2000) used Fleming et al's calculation of angular momen- 
tum transport as a function of Re m in a time-dependent 
numerical model for DN outbursts. He found that, as 
proposed by GM, MHD turbulence dies away in quies- 
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cence. This interplay between the thermal ionization in- 
stability and the magneto-rotational instability provides 
a physical motivation for the observational inference that 
c*coid <C ahot; it suggests that either accretion completely 
stops in quiescence or an additional, less efficient, angular 
momentum transport mechanism operates (perhaps spi- 
ral density waves excited by the companion; Spruit 1987; 
Menou 2000). 

Thin accretion disks around supermassive black holes 
in active galactic nuclei are also expected to have regions 
with T c ~ 10 4 K, and are thus also subject to the ther- 
mal ionization instability (e.g., Lin & Shields 1986; Can- 
nizzo 1992; Siemiginowska, Czerny, & Kostyunin 1996). 
Applications of the disk instability model to AGN have, 
however, received less attention than those to DN and 
SXTs; this is no doubt because of the much longer recur- 
rence times, which are observationally inaccessible in in- 
dividual sources. The thermal ionization instability could, 
however, significantly influence the luminosity function of 
AGN, perhaps accounting for the short "duty cycle" of 
bright Seyferts and quasars (e.g., Mineshige & Shields 
1990; Siemiginowska & Elvis 1997). In addition, it could 
lead to a radially varying accretion rate in the disk, com- 
plicating comparisons between theoretical models and ob- 
served broad band spectra. Understanding the role of the 
ionization instability and other sources of large amplitude 
variability in AGN is thus of considerable importance. 

In this paper we generalize GM's analysis and argue 
that, in contrast to SXTs and DN, the gas is well coupled 
to the magnetic field in AGN disks even on the cold branch 
of the thermal limit cycle (quiescence). This is essentially 
because the disks in AGN are physically much larger than 
those in DN; for a fixed resistivity, the magnetic Reynolds 
number is thus much larger. Our analysis implies that the 
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thermal ionization instability is unlikely to be a significant 
source of large amplitude variability in AGN. 

Accretion disks around supermassive black holes also 
differ from disks in DN and SXTs in that they may develop 
both local and global gravitational instabilities, which can 
lead to efficient angular momentum transport (e.g. Shlos- 
man, Begelman & Frank 1990; Kumar 1999; Gammie 
2000). In this paper, we argue that the degree to which 
the magneto-rotationally and gravitationally unstable por- 
tions of the accretion disk overlap is important for under- 
standing time-dependent accretion in AGN. 

We proceed as follows. We first calculate the structure 
of AGN disks in regions subject to the thermal ionization 
instability (§2). We then estimate the importance of re- 
sistivity, ambipolar diffusion, and Hall currents in these 
regions (§3). In §4 we reassess the role of the ionization 
instability in AGN. In §5 we discuss the interplay between 
MHD (un)stable 2 and gravitationally (un)stable regions of 
AGN accretion disks, and in §6 we conclude. 

2. THE LOCAL THERMAL STRUCTURE OF AGN DISKS 

We consider a gaseous, thin accretion disk around a su- 
permassive black hole (BH) of mass M. We ignore X-ray 
irradiation in calculating thermal equilibria (but see §5). 3 
The local thermal equilibria of the disk are then found by 
equating the local viscous dissipation rate Q + to the radia- 
tive cooling rate Q~ , where Q~ is a function of the disk 
opacity. Such equilibria have previously been calculated 
by several authors for AGN disks (see, e.g., Mineshige & 
Shields 1990; Cannizzo 1992; Siemiginowska et al. 1996). 
As for disks in DN and SXTs, these equilibria are ther- 
mally and viscously unstable in regions of partial ioniza- 
tion, corresponding to T c ~ 10 4 K. 

Figure |l| shows two thermal equilibrium curves ("S- 
curves") for an annulus located at R = 10 14 cm around 
a BH of mass M — 10 7 Af Q ; for the upper S-curve, 
(Jhot = a coid = 0.1, while for the lower S-curve, ahot = 
0.1 3> a co id = 0.01. The thermal equilibria have been 
found by calculating a grid of models for the disk verti- 
cal structure, for different values of the surface density £ 
and the central temperature T c , at the radius of interest 
(Hameury et al. 1998). Equilibria correspond to (S,T C ) 
pairs such that Q + = Q~. 

The middle, inverted sections of the S-curves are where 
the gas in the disk is partially ionized and unstable. Any 
disk annulus ostensibly in this state instead undergoes a 
limit cycle in which it time-dependently evolves between 
the two stable (hotter and colder) thermal equilibria. 

Point A in Figure 1 represents the lowest temperature 
reached by any unstable annulus during the thermal limit 
cycle. This point is of interest because it is the place in 
the limit cycle where the gas in the disk is least coupled to 
magnetic fields (it has the smallest ionization fraction in 
LTE). As a conservative estimate, we use T c « 2000 K for 
the central temperature of the unstable annulus at point 



A. This is appropriate for a co id = 0.01. 4 If, as we argue in 
§3-4, MHD turbulence is responsible for accretion in AGN 
disks in the vicinity of point A, a co id ~ ahot ~ 0.1 would 
be more appropriate (see Hawley 2000 for global MHD 
simulations which give a s=s 0.1). This only strengthens 
our conclusions. 

3. COUPLING OF THE GAS TO THE MAGNETIC FIELD 

GM showed that Re M ~ 3 x 10 3 at point A in DN disks; 
here we generalize their analysis and assess the coupling of 
gas to magnetic fields in unstable annuli of accretion disks 
around compact objects of arbitrary mass. 

The importance of ohmic diffusion can be estimated us- 
ing the magnetic Reynolds number: 
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where H 
the sound 

ity, Ven ~ 



a c s /Q is the scale height of the disk, c s is 
speed rj = c 2 m e ^ e „/(47rn e e 2 ) is the resistiv- 
8.3 x 10 _10 T 1 / 2 n„ s^ 1 is the electron-neutral 
collision frequency (Draine, Roberge, & Dalgarno 1983), 
and n n (n e ) is the number density of neutrals (electrons). 
For the relatively low density disks of interest in AGN, 
ambipolar diffusion can be more important than resistive 
diffusion. Its importance can be estimated using the di- 
mensionless number (e.g., Blaes & Balbus 1994) 
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ReA measures the ratio of the neutral-ion collision fre- 
quency, v n i, to the orbital frequency in the disk, fi; for 
ReA ^S> 1, the neutrals are well coupled to the magnetic 
field by collisions with the ions. In equation (||), v n i = jpi, 
7 w 3.5 x 10 13 cm 3 s" 1 g" 1 is the drag coefficient for ion- 
neutral collisions 5 (Draine et al. 1983), and pi is the ion 
mass density. 

To estimate ReM and ReA, we use the physical condi- 
tions in the limit cycle least favorable to MHD turbulence, 
namely those at point A; in particular, we use Ta ~ 2000 
K, a somewhat conservative value (§2). From our numer- 
ical calculations of S-curves, we find that the surface den- 
sity at point A can be approximately fit by (see also Can- 
nizzo 1992; Siemiginowska et al. 1996) 



T,a « S a 01 
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where £ ~ 1000 g cm~ 2 , ao.i = (a/0.1), r is the ra- 
dius of the unstable annulus in Schwarzschild units, and 
mg — (M I 10 9 Mq) . Finally, for a given rate at which mass 
is fed to the accretion flow at large radii (M, in Mq yr _1 ), 
the radius around which the ionization instability sets in 
is given by 

-2/3 



TA 



r M 0A rr,„ 



(4) 



2 Throughout this paper, we use "MHD unstable" to mean unstable to developing MHD turbulence via the magneto-rotational instability. 

-ray irradiation can lead to higher temperatures than considered below; this can substantially change the thermal ionization instability 
(e.g., Dubus et al. 1999). 

This result can be checked by solving the vertically-averag ed equation Q+ = (9/8)i^Ef2 2 = Q~ = 4o"T c 4 /(3Sk) for an optically-thick disk, 
where v is the kinematic viscosity and k is the opacity (which can be obtained from a standard opacity table for the range of density and 
temperature of interest). It is also interesting to note that the temperature Ta — 2000 K at point A on the S-curve is, for a co id ~ 0.01, a 
generic prediction for any black hole mass and at any radius in the disk (see, e.g., GM for an identical result in the case of a DN disk). 

5 We take rrii = 30mu and m n = 2.33mj^, where rrii is the ion mass, m n is the neutral mass, and mu is the mass of the hydrogen atom. 
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where r « 300. Equations (||) and (Q) can be roughly 
reproduced analytically by fixing T ~ 10 4 K in a Shakura- 
Sunyaev disk (thus isolating the portions of the disk where 
hydrogen starts being partially ionized). It is also inter- 
esting to note that equations (^) and (||) are applicable 
to DN and SXTs as well as AGN, with m 9 then ~ 10~ 9 , 
M ~ 10~ 9 M Q yr -1 , and r A ~ 10 5 . 

Using the Saha equation for a solar metallicity gas, we 
find that, in LTE, the ionization fraction at Ta ~ 2000 K 
in the thermally unstable annulus is x e = n e /n n « 
10~ 6 - 10~ 7 for M e [1, 10~ 10 ] M yr- 1 . 6 Substituting 
equations (||) and (||) into equations (|l|) and @ we then 
find the smallest values of ReM and ReA in the thermal 
limit cycle (obtained at point A): 



Re M ~ 10 ! 



1 300/ 



3/2 



and 



Re A w 3 x 



10 7 a, 



3/4 
0.1 



V 300/ 



1000 / 



(5) 



(6) 



For a given accretion rate, equations (||) and (||) determine 
the coupling of the gas to the magnetic field in the annulus 
which is unstable to the ionization instability; this annulus 
lies at a radius given by equation (0). Note that Reu and 
ReA are independent of black hole mass at fixed M . 

Hall currents can significantly modify the dynamics of 
low density, weakly-ionized disks (e.g., Wardle 1999; Bal- 
bus & Terqucm 2000). Their importance can be measured 
using the ion Hall parameter, (3i = uji/"fp n , where u>i is 
the ion cyclotron frequency. For /3, 3> 1 the ambipolar 
diffusion limit applies, while for /3i <§; 1 Hall currents are 
important. Assuming equipartition magnetic fields, we es- 
timate A « 0.3M 01 (S /1000)- 1/2 in the unstable annuli 
of interest. The role of Hall currents in modifying MHD 
turbulence is determined by Ren ~ PiReA, since Pi ~ 1, 
ReA is a reasonable proxy for Ren- if ReA S> 1, then 
Ren 3> 1, MHD turbulence is fully developed, and Hall 
currents are unimportant. By contrast, if ReA ~ 1 MHD 
turbulence is modified (and Hall currents must be consid- 
ered because of their potentially destabilizing effects). 

4. IMPLICATIONS 

MHD turbulence and its associated angular momentum 
transport are expected to be suppressed for ReM and/or 



By contrast, for reasonable M in AGN, ReM 3> Re M 
and Rca 3> Re c A . This implies that the disk is MHD- 
turbulcnt throughout the thermal limit cycle. In con- 
trast to DN and SXTs, the efficiency of angular momen- 
tum transport should then be comparable on the hot and 
cold branches (a co id ~ ai 10 t)- In this case global time- 
dependent calculations show that the thermal ionization 
instability is probably not a significant source of time- 
dependent accretion in AGN; it leads to luminosity "flick- 
ering," rather than large amplitude outbursts analogous 
to DN and SXTs (e.g., Smak 1984; Mineshige & Shields 
1990). 

The analysis presented above does not, of course, pre- 
clude that accretion in AGN is highly time dependent; it 
does, however, suggest that the thermal ionization insta- 
bility, as applied to DN and SXTs, is not relevant to AGN. 
Radiation-pressure induced thermal and viscous instabil- 
ities at small radii may be an important source of time 
dependence in AGN accretion disks (Lightman & Eardley 
1974; Piran 1978). Another, less well explored, possibility 
lies in the interplay between MHD (un)stable and grav- 
itationally (un)stable portions of the accretion disk. We 
briefly address several of the relevant issues below; much 
more theoretical work on these problems is needed. 



5. MAGNETO-ROTATIONALLY AND GRAVITATIONALLY 
UNSTABLE REGIONS 

5.1. Magneto-rotational Stability 

An obvious shortcoming of the argument in §3 is that 
it is local. We examined annuli in thin accretion disks 
which are subject to the thermal ionization instability and 
found that they are always MHD turbulent so long as 
M 3> 10~ 7 M Q yr _1 . At radii somewhat larger than the 
thermally unstable region, however, the disk is everywhere 
stable, but on the cold branch of the S-curve in Figure 1; if 
viscous dissipation is the only source of heating, the cen- 
tral temperature in the disk eventually drops to <C 1000 
K. In LTE the ionization fraction then becomes very small 
and MHD turbulence dies away. 

X-rays and cosmic-rays may provide sufficient heating 
and/or nonthermal ionization to maintain magnetic cou- 
pling in the disk at large radii. Blaes & Balbus (1994) 
estimated ReM 3> 1 and ReA ~ 1 from cosmic-ray ioniza- 
tion in the parsec scale circumnuclear disk in our Galaxy. 
Simple estimates of X-ray ionization yield similar results. 

As a concrete example, consider the w 0.1 — 0.25 pc mas- 
ing disk around the M = 3.6 x 10 7 M Q black hole in NGC 
4258 (Miyoshi et al. 1995). Neufeld & Maloney (1995; 
NM) proposed that the masing region is strongly irradi- 
ated by the central X-ray source, due to a significant warp 
in the disk. For a disk illuminated at an angle cos -1 /! to 
the normal direction, the incident X-ray flux on the disk 
is / w (^L X )I{^R 2 ), where L x ~ 4 x 10 40 ergs s" 1 for 
NGC 4258 (Makishima et al. 1994) and fi may be as large 
as « 0.2 (NM; Herrnstein 1997). The volume ionization 
rate in the bulk of the disk due to the incident X-ray flux 
is f w mm(l,T)f/(HE l ) cm" 3 s" 1 , where J5j w 37 eV 
is the typical energy required to generate one secondary 
electron (Glassgold, Najita, & Igea 1997), r = /i _1 S/S c is 

6 The neutral density at point A in the unstable annulus is n n ~3x 10 15 M — 0,2 cm -3 . From the Saha equation, x e oc n n 1,/2 oc M 0,1 . Thus, 
at fixed T c , the 10 order of magnitude change in M from DN to luminous AGN implies a change in x e by a factor of f» 10. Larger values of 
x e are expected in high metallicity AGN disks because n e is dominated by contributions from metals with small ionization potentials. 



ReA less than critical values = Re 



M 



and Re c A . 



Numeri- 



cal simulations with explicit resistivity are consistent with 



Re 



M 



2000), 
a 100 
1995; 



10 4 (Hawley et al. 1996; Fleming et al. 
while simulations of ion-neutral disks give Re c A 
(Hawley & Stone 1998; see also Mac Low et al. 
Brandenburg et al. 1995). 

Equations (§) and (||) then show that there is an im- 
portant difference between accretion in DN and SXTs, 
where M ~ 10" 9 Moyr -1 , and accretion in luminous 
AGN, where M ~ 0.01 - 10 M yr -1 . As shown by GM, 
Re M £ Re c M in DN and SXTs. MHD turbulence is thus 
suppressed on the "cold" branch of the thermal limit cycle. 
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the shielding optical depth of the disk, and E c w 1 g cm 
is the column density corresponding to order unity optical 
depth for hard X-rays ib 10 keV (the factor min(l, t) is in- 
cluded to take into account the possibility of a disk which 
is optically-thin to X-rays). NM find n„ ps 3 x 10 7 cm~ 3 
in the masing region and so £ ps 0.1 g cm -2 . Balancing 
X-ray ionization against dissociative recombination at a 
rate 8.7 x 10- 6 n 2 T~ 1/2 cm" 3 s" 1 yields n e w 300 cm" 3 
and thus x e ~ 10~ 5 . From equations (Q) and (||), this 
corresponds to Reu ~ 10 11 and i?eA ~ 10 3 . Somewhat 
counterintuitively, the disk in NGC 4258 may be both mas- 
ing and MHD-turbulent! Observational support for this 
conclusion is provided by the work of Wallin, Watson, & 
Wyld (1998, 1999); they showed that the maser spectral 
line features observed in NGC 4258 can be reproduced 
reasonably well by including a turbulent velocity field in a 
thin Keplerian disk model. 

In order to account for the absence of masing at radii 
< 0.1 pc, NM propose that X-ray irradiation "shuts off" 
inside ps 0.1 pc (because the warp decreases with radius, 
/U — > 0). If correct, MHD turbulence probably ceases 
to be a viable angular momentum transport mechanism 

(cosmic-ray ionization alone yields Rem ~ lO 7 ^ 2 , but 

ReA ~ O.l^ 2 , where £ = 10~ 17 £i7 s _1 is the cosmic-ray 
interaction rate per hydrogen atom). 7 

This illustrates a more general worry; the role of non- 
thermal heating and ionization is quite uncertain, depend- 
ing on, e.g., the geometry of the disk (warps, flaring, etc.), 
the presence of X-ray scatterers, and the possibility of an 
anomalously large cosmic-ray flux due to a nearby jet. It 
is possible that many parsec scale disks in AGN are not 
magnetically well coupled. In this case gravitational insta- 
bility is the most promising angular momentum transport 
mechanism. 

5.2. Gravitational Stability 

Thin disks are locally gravitationally unstable for Q <, 1, 
where Q = c s Sl/ (hGYj). The nonlinear outcome of local 
gravitational instability depends on the ratio of the cool- 
ing time to the rotational period of the disk (Shlosman & 
Begelman 1989; Shlosman et al. 1990; Gammie 2000). For 
long cooling times, the disk reaches a steady state with effi- 
cient angular momentum transport produced by "gravito- 
turbulence" (Shlosman et al. 1990; Gammie 2000). For 
parsec-scale disks in AGN, however, the cooling time is 
generally much shorter than the rotational period of the 
disk; in this case the disk fragments, and probably under- 
goes efficient star formation (Shlosman & Begelman 1989; 
Gammie 2000). It is unlikely that local gravitational in- 
stability produces significant accretion in this regime. 

For Mdisk ~ 2irYjR 2 > M, however, the disk is glob- 
ally gravitationally unstable. Using M — 2irY,R\v r \ and 
| TV | ~ ac 2 /«0, where v<f, is the (Keplerian) rotational ve- 
locity, this condition can be rewritten as (Shlosman & 
Begelman 1989) 

2 

M J> ~ KTViTaWioo M Q yr" 1 , (7) 

where T = 10 3 T 3 K and vj, = 100wi O o km s" 1 . 

7 The ion Hall parameter is ft ps 10 3 assuming equipartition fields, 
vicinity of the AGN could, however, provide significant ionizing phot' 
8 Accretion may proceed in an ionized surface layer analogous to C 



If nonthermal heating and ionization are insufficient to 
maintain MHD coupling, disks will become MHD stable 
if T 3 <> 1. For luminous AGN with M £ 10" 2 A^yr" 1 , 
however, the disk is globally gravitationally unstable in 
this region (eq. Q). It is therefore likely that bar forma- 
tion in the gaseous disk leads to significant angular mo- 
mentum transport - and thus accretion - on a dynami- 
cal timescale (Shlosman et al. 1990). At smaller radii, 
where T3 J> 1, the disk is both MHD and gravitationally 
unstable; the dynamics in this regime has been poorly ex- 
plored but may be accessible with numerical simulations. 
More generally, in luminous AGN the accretion disk proba- 
bly evolves from gravitationally unstable and MHD stable 
at large radii to gravitationally unstable, MHD unstable, 
and finally to gravitationally stable, MHD unstable at the 
smallest radii. It is rather uncertain whether these transi- 
tions occur at constant M. 

In low- luminosity AGN with M <, 10~ 2 M Q yr~ x equa- 
tion (Q) shows that the disk is globally gravitationally sta- 
ble where it is MHD stable; for M < 10" 3 Afgyr" 1 the 
disk has Q J> 1 and is locally stable as well; this raises 
the interesting possibility of an MHD and gravitationally 
stable "dead zone" in which accretion may proceed very 
inefficiently (or perhaps not at all). 8 Mass will build up in 
this region until Mdisk ^ M, until the disk warp evolves so 
as to provide significant X-ray ionization, or until a global 
hydrodynamic instability sets in (analogous to GM's pro- 
posal for DN). It is possible that such piling up of mass is 
presently occurring in many nearby galaxies. In particular, 
returning again to NGC 4258, NM's masing model gives 
M ps 10~ 5 ao.i M yr _1 and thus a very gravitationally 
stable disk: M disk ~ lOOAi© < M and Q « 100 > 1 (see 
Gammie, Narayan, & Blandford 1999 for alternate M es- 
timates) . If X-ray irradiation and MHD turbulence indeed 
shut off for R < 0.1 pc, it is unlikely that significant ac- 
cretion occurs at these radii. It is then also unlikely that 
NM's accretion rate estimate is relevant for the energy- 
producing portions of the accretion flow at small radii. 

6. CONCLUSION 

We have studied the combined role of the thermal ion- 
ization, magneto-rotational and gravitational instabilities 
in thin accretion disks around supermassive black holes. 

In contrast to dwarf novae and soft X-ray transients, we 
find that there is no physical motivation for a reduced ef- 
ficiency of angular momentum transport when AGN disks 
become neutral via the thermal ionization instability. This 
is because even the predominantly neutral gas is well cou- 
pled to the magnetic field (Re\i 3> 1 and Re a 1). 
MHD turbulence due to the magneto-rotational instability 
should thus be present and similar to that in an ionized 
plasma. In this case, numerical simulations have shown 
that the thermal ionization instability leads to small am- 
plitude luminosity "flickering," rather than large ampli- 
tude outbursts analogous to DN and SXTs. 

We propose, however, that accretion in AGN can be 
time-dependent because of the interplay between the 
magneto-rotational and gravitational instabilities, both of 
which are viable candidates for transporting angular mo- 
Hall currents should not change this result. Scattered X-rays from the 
; their importance is very uncertain, 
imie's (1996) proposal for T-Tauri disks. 
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mcntum in AGN accretion disks. The hot inner regions 
of the disk are preferentially subject to the magneto- 
rotational instability, while the cold outer regions are 
prone to developing gravitational instabilities. In high lu- 
minosity systems (accretion rate M <; 10~ 2 M Q yr -1 ), the 
magneto-rotationally and gravitationally unstable regions 
of the accretion disk overlap, while in lower luminosity 
systems (M <, 10~ 3 M yr _1 ) they may not. In systems 
where these regions do not overlap, such as nearby low- 
luminosity galactic nuclei on scales of ~ 0.1 — 1 parsec, 
there is no well understood process for removing angular 
momentum from accreting gas. Mass may therefore be 



temporarily building up in a non-accreting "dead-zone." 
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INSTABILITIES IN SUPERMASSIVE BLACK HOLE DISKS 
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Fig. 1. — Example "S-curves" (curves of thermal equilibria in a surface density - central temperature (S-T c ) diagram) for a thin accretion 
disk around a supcrmassive black hole. The disk annulus is located at R = 10 14 cm around a black hole of mass M = 10 7 Mq. For the upper 
S-curve, the values of the viscosity parameter are a^ ot = a co id = 0.1, while for the lower S-curve a^ ot = 0.1 2> Q co id = 0.01. The point 
labeled A on each of the lower branches corresponds to the coldest state reached by the annulus during its thermal limit cycle evolution; this 
is where the gas is least coupled to the magnetic field. 



